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ABSTRACT 

Combining measurements of the galaxy power spectrum and the cosmic microwave 
background (CMB) is a powerful means of constraining the summed mass of neutrino 
species ^m„, but is subject to systematic uncertainties due to non-linear structure 
formation, redshift-space distortions and galaxy bias. We empirically test the robust- 
ness of neutrino mass results to these effects by separately analyzing power spectra 
of red and blue galaxies from the Sloan Digital Sky Survey (SDSS-II) Data Release 
7 (DR7), combined with the CMB five-year Wilkinson Microwave Anisotropy Probe 
(WMAP5) data. We consider fitting for a range of maximum wavenumber k using 
twelve different galaxy bias models. For example, usi ng a new model b ased on pertur- 
bation theory and including redshift space distortions (jSaito et al. 2009), the all-galaxy 
power spectrum combined with WMAP5 for a wavenumber range of k < 0.2/iMpc -1 
yields 95% CL ^to„ < 0.46 eV. The red and blue galaxy power spectra give 0.41 
and 0.63 eV respectively for this model. Using mock catalogues, we find the expected 
difference in these limits assuming a true neutrino mass of zero is 0.10 ± 0.14 eV. 
Thus the difference of 0.22 eV between upper limits on neutrino mass for red and 
blue galaxies is approximately la from the expected value. We find similar results for 
the other models and k ranges tested. This indicates good agreement for current data 
but hints at possible issues for next-generation surveys. Being able to perform such 
systematic tests is advantageous, and future surveys would benefit by including broad 
galaxy populations and luminosities that enable such a decomposition. 

Key words: neutrinos - cosmological parameters - cosmology: observations - galax- 
ies: statistics - cosmic background radiation - large-scale structure of Universe. 



1 INTRODUCTION 

Atmospheric and solar neutrino experiments show that neu- 
trinos have mass and that there is significant mi xing between 



the various neutrino interaction eigens t ates (iDavis 



2005 



119981: I Ahmad et all [20021 : lAshie et al 
IWendell et ajJ|201C| ). Cosmological measure- 



1996 



Fukuda et al. 
Martini |20ol 

ments nicely complement these experiments by constraining 
the sum of the masses of the different eigenstates m Vi . 
The cosmological experiments exploit differences between 
the behaviour of cold dark matter (CDM) and massive neu- 
trinos, which free-stream when relativistic, removing small- 
scale structure that would have for med in a CDM-only uni- 
verse (for recent reviews see, e.g., lElgar0v fc LahavlBuOl 
iLeseourgues fc Pastor] 120061 : iHannestad 1201(f ). To first or- 
der, the suppression depends on Q. v , the total mass density 
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comprised of neutrinos relative to the critical density of the 
Universe. This is related to the sum of the mass eigenstates 
by X] m v = (94.1 fo 2 eV), where h is the Hubble constant 
Ho divided by 100 kms _1 Mpc _1 . Thus cosmological mea- 
surements primarily probe the total particle mass summed 
over neutrino species. In principle, the neutrino mass hierar- 
chy could be measured by resolving the free-streaming scale 
of individual neutrino species, but curren t experiments do 
not have the sensitivity required to do th is |Hannestadll2003l : 
iTakada et"ai1l2006l ; [jimenez et alJl2010T ). 

The most recent upper limit on the summed neutrino 
mass published by the Wilkinson Microwave Anisotropy 
Probe (WMAP) tea m is £m„ < 0.58 eV at 95% CL 
(|Komatsu et aljhoiti ) .This limit combines CMB data from 
the 7 year WMAP data release (WMAP7) with a measure- 
ment of Ho from type la supernovae (SN la) and the Baryon 
Acoustic Oscillation (BAO) information from the 2-degree 
Field Galaxy Redshift Survey (2dFGRS) and Sloan Digital 
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Sky Survey (SDSS-II). Adding large scale structure infor- 
mation from surveys of galaxies, galaxy clusters, or weak 
lensing pushes this limit down even further: other groups 
have combined WMAP5 with SN la, BAO, and large scale 
structure information a nd derive limits as low as Y] m v < 
0.28 eV at 95% CL (e.g- JVikhlinin et al.ll2009| : |Tereno et al 



2009; iReid et all l2010al ; iMantz et al.l I2OI0I ; iThomas et al 



2010J). The challenge now is to reliably bring down the up- 



per limits to the 0.1 eV level or even detect the neutrino 



data (|Lesgourgues et al. 


2004; Lesgourgues & Pastor! l2006t 


lHannestad & Wondl2007 


:lLahav et alJboioh. 



Unfortunately, our ability to measure the summed neu- 
trino mass using large scale structure information such as the 
galaxy power spectrum will likely be limited by our lack of 
knowledge about other effects that change the expected am- 
plitude of galaxy clustering on small-scales. These effects in- 
clude the formation of non-linear matter structures and the 
non-linear galaxy peculiar velocities within these structures. 
The peculiar velocities affect the measured power spectrum 
because we infer distances from redshifts, and peculiar ve- 
locities are misinterpreted as being due to the Hubble flow. 
Thus the random motion of galaxies within collapsed struc- 
tures acts to damp the power spectrum on small scales. 

In addition, galaxies are not expected to exactly Pois- 
son sample the distribution of matter in the Un iverse 
l|Sheth fc Lemsonl Il999l ; [c asas-Miranda et al. I 120021 ); evi- 



dence for this has been seen in, e.g., IWake et al. {200I). 
This severely complicates the extraction of cosmological in- 
formation from galaxy surveys. Furthermore, it has been 
known for some time that different populat ions of galaxies 
demonstrate differen t clustering strengths (|Davis & Geller 
19761 : iDresslerl [l980l; iPark et al.1 1 19941; iPeacock fc Doddsl 



1994; ISeaborne et al.Ml999l; iNorberg et all l200ll. |2002|; 



Zehavi et al .1120021. 12005; Wild et al.ll2005; Collist er fc Lahavl 



20051 ; iLi et al. 2006; Swa nson et~ai1l2008h . showing that they 



cannot all have a simple relationship linking their distribu- 
tion with that of the matter. This galaxy bias (the rela- 
tionship between the galaxies and the matter) is the most 
pernicious physical effect that limits o ur ability to extract 
cosmological data from galaxy surveys (jPercival et al ] |2007l ; 
ISanchez fc Coiell200Sl ). 

Because of the uncertainty regarding these effects, many 
cosmologists (including the WMAP team) choose not to use 
the galaxy power spe ctrum in making their cosmological 
param eter constraints (iDunklev et al . 2009 1; iKomatsu et alj 
120101 ) and instead encapsulate the information from galaxy 
surveys as a prior on the geometrical distance to the mean 
redshift of the galaxy survey based on using the BAO acous- 
tic peak as a standard ruler. This is effective but discards 
a great deal of information contained in the shape of the 
power spectrum that is especially valuable for making neu- 
trino mass constraints. In order to fully exploit the informa- 
tion contained in the galaxy distribution, we must develop 
a more sophisticated understanding of how the power spec- 
trum is affected by galaxy bias, non-linear evolution, and 
redshift space distortions. Of these three systematic effects, 
the galaxy bias is in some sense the most troublesome: non- 
linear evolution and redshift space distortions can be studied 
in N-body simulations allowing only for evolution through 
gravity, but modelling galaxy bias in detail requires hydro- 



dynamic simulations or semi-analytic galaxy formation mod- 
els. 

The aim of this paper is to quantify the impact of these 
systematic effects on measurements of the summed neutrino 
mass from the galaxy power spectrum, focusing especially on 
the effects of galaxy bias. We use data from the final data re- 
lease of the SDSS-II survey, which we split into the bimodal 
populations of blue and red galaxies. By comparing results 
from two galaxy populations with differing properties, we 
can study the impact of astrophysical effects while bypass- 
ing the need to model biasing in detail. We also consider 
how the range of wavenumber k fitted affects the results. 
Using data at larger k values adds more information about 
neutrinos from the effects on the small-scale structure, but it 
comes at the price of requiring more complicated theoretical 
modelling. We fit our subsamples and the combined sample 
with 12 different power spectrum models and combine with 
CMB constraints from WMAP5. 

This work is timely for two reasons. Firstly, the up- 
coming generation of galaxy surveys s uch as the Dark En- 
ergy Survey (DES; lAbbott et al.ll2006h . the Panoramic Sur- 
vey Telescope fc Rapid Response System (Pan- STARRS; 
iMorgan fc Kaiserll2008l ; iMorgan et apbOOSl'), and the L arge 
Synoptic Survey Telescope fLSST; lAbell et all 120091 ) will 
significantly reduce the statistical uncertainty on power 
spectrum measurements, making a detailed understanding 
of the systematic uncertainties increasingly essential. Sec- 
ondly, we can now make use of new perturbation theory re- 
sults that provide more accurate models for th e non-linear 
matter power spectrum including neutrinos dSaito et al.l 
l2009l ; lTaruva et al.|[2009l ): along with lSaito et all (|2010h this 
analysis is one of the first applications of these models to ob- 
servational data. 

This paper is organised as follows. We start by describ- 
ing the data that we use in Section [5] The method by which 
these data are analysed, including power spectrum calcula- 
tion and likelihood analysis is described in Section[3] Results 
are presented in Section fj] and are discussed in Section [S] 
For simplicity, here we assume a cosmological model of flat 
ACDM plus massive neutrinos, but we note that alternative 
cosmological models (in particular those with dark energy 
equation of state w 7^ — 1) can impact neutrino mass mea- 
surements as wel l (|Hannestadl 120051 ; iKiakotou et alj 120081 ; 
IReid et al.lboiOal 'l and will be investigated further in future 
work. Specifically, we fix Q,m = 1, the dark energy equa- 
tion of state w — —1, the tensor amplitude and running of 
the spectral index to zero, and the effective number of rel- 
ativistic degrees of freedom iV c ff = 3.04. We convert galaxy 
redshifts to distances using Q,m = 0.25 and Qa = 0.75, and 
assume that this is the true cosmology when estimating er- 
rors on these measurements. 



2 DATA 

2.1 Galaxy data 

The SDSS-II experiment used a 2.5m telescope l|Gunn et alj 
2006), to obtain a spectroscopic sample of galaxies selected 
to a limiting Galactic extinction-corrected Petrosian mag- 
nitude r < 17.77, or r < 17.5 in a small su bset of the 
early data from the survey (jStrauss et all [2002). There are 
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Figure 1. Histogram of the probability distribu tion for Y*, m v 
reproduced using the WMAP5 MCMC chain from lDunklev et al.l 
11200911 . 



approximately 90 galaxies per square de gree, with a me- 
dian redshift of z = 0.1. The DR7 sample (|Abazaiian et al] 
l2009h used in our a nalysis includes 669 905 main galaxies 
l|Strauss et al.ll2002h . We correct for fibre packing problems 
where close galaxies cannot both be observed by assigning 
the redshift of the nearest observed galaxy to a galaxy which 
was no t observed, matching the methodology of Zehavi et al.l 
(2002). This is the sub sample of SDSS main galaxies used 
bv lPercival et all (|2010T l. Galaxy redshifts were converted to 
distances using our fiducial cosmology (flat ACDM model 
with tt m = 0.25). 

We have fc-correcte d galaxy lum i nosities using the 
methodology outlined in iBlanton et al.l (|2003al lbh . We also 
use a z = 0.1 shifted r-band filter to defi ne our lumi- 
nosities (as discussed in IBlanton et al] 2003a, b), which we 
refer to as Mo.i,.. Absolute magnitudes and fc-corrections 
were calculated assuming Ho = 100kms _1 Mpc~ 1 , and our 
fiducial cosmology, and we have applied the recommended 
AB corrections to the observed SDSS magnitude system 
l|Smith et al.ll2002l ). 

This sample of galaxies was split into red and blue sub- 
samples defined by a constant colour cut of Mo.i g — Mo.i r = 
0.8. This approximately divides the bi modal distribution of 
galax ies into red and blue types (e.g. ICresswell fc Percivall 
2009). Our results are not sensitive to the exact nature of 
this division. The red and blue subsamples have median red- 
shifts of z = 0.111 and z = 0.085 respectively. 



ning spectral index, or dark energy equation of state w. As 
discussed in llchikawa et alj (|2005h , the best one can hope 
to do for the neutrino mass limit using only CMB data is 
^m„ < 1.2 eV. Since the WMAP5 limit is already nearly 
at this value, we do not expect significant improvement by 
going to WMAP7. This is borne out by t he fact that the 
WMA P7-only limit is also J2 m " < 1-3 eV jKomatsu et all 
l20ld) . 

We have downloaded th e WMA P5 ACDM+M„ MCMC 
chain from iDunklev et"aH (|2009l ) from the LAMBDA 
archival} This model uses eight parameters to fit to the 
CMB temperature and polarization power spectra: the cold 
dark matter density Q c h 2 , the baryon density Q.bh 2 , the neu- 
trino mass density Q^h 2 , the dark energy density Qa, the 
amplitude of curvature fluctuations A|>, the scalar spectral 
index n s , the reionization optical de pth t, and the ampli- 
tude Asz of a Sunyaev-Zel'dovich (SZ. ISunvaev fc Zeldovichl 
1 1970T ) fluctuation template spectum. We express these pa- 
rameters in a more digestible form using the flatness as- 
sumption, ~^2,m v = Q v (94.1 h 2 eV), and the definition of 
the linear theory density fluctuation amplitude at 8 ft" 1 Mpc 
scales: 



0"8 



2^2/ [3 (sins- 



x cos x) /x 3 ] 2 Pii n (k) k 2 dk 



1/2 



(1) 



with x = k x 8/i _1 Mpc. This gives a parameter set l = 
(yih, Qb/Qm, h, as, Us, X] m ") P ms the CMB-only parame- 
ters (t, Asz) with f2 m = Q c + Qt + fl v . 

Using this chain we have reproduced the ^ m u < 1.3 eV 
limit, illustrated in Fig. [1] This is a histogram of the N = 
7503 values of "^2 m v = Mv in the MCMC chain weighted by 
the weights uicmb given in the chain. The 95% confidence 
limit on ^ m„ , which we denote as Ml lul , is defined in the 
usual way as a region enclosing 95% of the histogram area: 



.95 = 



J^M'KM}}™ W CMB 
L-ii=\ "'CMB 



(2) 



where i is an index running through the entries in the 
MCMC chain. 



3 ANALYSIS METHOD 

We explore the robustness of a neutrino mass limit produced 
by combining the CMB with information from the galaxy 
power spectrum. We do this by performing the analysis us- 
ing power spectra of different subsamples of galaxies, namely 
red and blue galaxies, and determining whether the result- 
ing limits are c o nsiste nt with each other, in the spirit of 
ISanchez fc Cold (|2008l ). We also test the effect of changing 
the range of k values used from the galaxy power spectra and 
using different prescriptions for the theoretical modeling of 
the nonlinear power spectra. 



2.2 CMB data 

We use the WMAP5 CMB data as our b aseline for neu- 
trino mass constraints (|Dunklev et al.ll2009l ). The WMAP5 
CMB data alone provide a limit of ^ra v < 1.3 eV (95% 
CL), robust to within 10% to varying tensor amplitude, run- 



3.1 Calculating power spectra and uncertainties 

Power spectra were calculated for the red, blue and com- 
bined catalogs using the Fourier method of iFeldman et all 



1 http://lambda.gsfc.nasa.gov 
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111994 ). as applied by IPercival et alj (|2010h . In this method 
a weighted galaxy over-density field is defined and Fourier 
transformed, then the spherically averaged power is mea- 
sured. The shot noise due to the discreteness of the galaxy 
locations is subtracted, assuming the galaxies are a Poisson 
sampling of the density distribution. We use the luminos- 
ity (but not co l our) d ependent galaxy weights advocated by 
IPercival et al.1 (|2004l ). These will at least partially correct 
for power spectrum shape changes caused by variations in 
the galaxy bias over the redshift range of the survey due to 
flux-limited selection and galaxy evolution. 

The power spectrum for the combined catalog is not 
independent of the red and blue catalogs since it includes 
the same galaxies. The all-galaxy power spectrum P a n (fc) is 
related to the red and blue galaxy power spectra P rr (fc) and 
Pbb (k) by 

Pan (A;) = ftPn (k) + 2f r fbP lb (k) + / b 2 P bb (k) 

where f r and /j are the fractions of red and blue galaxies in 
the sample (f r + fb = 1) and P r b (k) is the cross power spec- 
trum between red and blue galaxies. Thus P a n (k) contains 
the information from the two subsample power spectra plus 
the cross power spectrum. 

Because of the survey mask (both angular and ra- 
dial), the measured powe r spectrum is a convolution of 
the t r ue power spectrum l|Hauser fc Peebles! Ii~973l ; iPeeblej 
Il980l ; IPercival et al.l l200"lh . A discretized window function 
w as used to quic k ly per form this convolution, as described 
in IPercival et al.l |2007h . In order to calculate the covari- 
ances between the power spectrum band-powe rs, we have 
creat e d 10 000 Log-Nor mal (LN) density fields (|Coles et all 
1 19991 ; I Cole et al. 2005) with a power spectrum for a flat 
ACDM cosmology with 6> l =(fi A = 0.7, Q. b /Q. m = 0.15, 
h = 0.7, n s — 1, = 0), and normalisation matched to 

the amplitude of the measured power. 

From each LN density field 8 (x), we draw red, blue, and 
combined galaxy catalogues as Poisson samplings matched 
to the galaxy bias and spatial distribution of the SDSS-II 
samples. That is, the mean A (x) of the Poisson distribution 
is given by X(x) — n (x) (1 + b8 (x)), where n{x) is the 
mean density including the effects of the survey mask and b 
is the linear galaxy bias at large scales (k < 0.1 ftMpc -1 ) for 
the red and blue galaxy samples. Catalogues were calculated 
on a (512) 3 grid with box length 4000 h' 1 Mpc. Unlike N- 
body simulations, these mock catalogues do not model the 
growth of structure, but instead return a density field with a 
log-normal distribution, similar to that seen in the real data. 
The mock catalogues also do not include redshift space dis- 
tortions. Power spectra for the mock catalogues were calcu- 
lated following the same method used for the SDSS-II data 
and were used to compute covariance matrices for each sub- 
sample. 

A subset of 200 of the mock catalogue power spectra 
was also used to estimate the expected sensitivity of the 
neutrino mass limit. We do this by applying our analysis 
described in Section r3.4l to the mock power spectra and de- 
fine a range using the mean ± one standard deviation of the 
95% CL limits on ^ m v from the mocks. Because the mock 
catalogues have ^ m v — and linear bias, the expected 
sensitivity range is the range in which we would expect to 
measure the neutrino mass limit assuming that the neutrino 
mass were equal to zero and that we understood bias per- 



fectly. This provides a baseline for interpreting the neutrino 
mass limits determined from the actual data. 



3.2 Models for galaxy bias and nonlinear 
evolution 

Our starting point for theoretical modeling is the 
linear power spectrum Pn n (k, z) calculated us- 
ing CAMrJf] for our set of cosmological parameters 
% = (Qa, r2(,/n m , h, as, n s , ^m v ). We evaluate P in (fc, z) 
at the median redshift of each galaxy sample (z — 0.111 for 
red, z = 0.085 for blue, and z = 0.099 for all galaxies) and 
drop the explicit dependence on 2 in our power spectra for 
notational convenience. 

We consider twelve different models for converting 
Pi n (k, z) into a galaxy power spectrum. The first six mod- 
els represent the typical models used in previous cosmolog- 
ical analyses of the neutrino mass and have either one or 
two additional free parameters, b and a free parameter that 
controls the amplitude of the non-linear change in shape be- 
tween galaxy and linear matter power spectrum, which can 
take on different values for different galaxy subsamples. 

For our remaining six models, we explore variations 
on new models for the matter power spectrum including a 
contribution from massive neutr i nos, based o n perturbation 
theory described in ISaito et al] i|200Sl , 120091 ); iTaruva et ail 
(2009). For this theory, the baseline model is 

P pt (fc) = /cbAtoop (k) + 2f Gh f v Pt v (k) + flPy^ (k) , (3) 

where / c b and f v are the mass fractions of cold dark mat- 
ter plus baryons and neutrinos respectively, P^ n (fc) is the 
linear theory power spectrum for neutrinos only, (fc) 
is the linear theory cross power spectrum between cold 
dark matter plus baryons and neutrinos, and Piio 0p (fc) is 
the nonlinear power spectrum of cold dark matter plus 
baryons computed to one loop in perturbation theory, 
which has been studied extensively in the literature (e.g., 
Juszkiewica Il98ll; iMakino et alj|l992; IJain k. Ber tschinger 
Il994l ; IScoccimarro fc Friemanl Il996l ; Heavens et all 11998 1 
Here we use 



Pitop (fc) = (fc) + P(2 b 2 ) (fc) + Pets) (fc) ■ 



(4) 



where P(22) (k) an d P(is) (k) are the approximations to the 
on e-loop corre c tion s pectra given by equations (32) and (36) 
in lSaito et ail l|2009h . If we replace Pn oop (fc) by Pfc (fc) in 
equation {(3j , the expression for P pt (fc) reduces to Piin(fc). 
The motivation for this model is that while cold dark mat- 
ter and baryon perturbations start to behave non-linearly at 
fc > 0.1 /iMpc -1 , perturbations of the neutrino density will 
tend to remain mostly linear due to the large velocity dis- 
pe rsion of the neutrinos. S i milar models have b een explored 
bv lHannestad et al.1 (120061 ) ; llchiki et al.1 (120091) and further 
improved bv lWond (|2008l ); iLesgourgues et al.l (|2009l ). 

The twelve models used in this analysis are as follows: 

(i) The simplest model for the galaxy power spectrum is 
that of constant linear bias of this linear power spectrum: 



Pgal (fc) = b 2 P in (fc) : 



(5) 



http : //camb/inf 0/ 
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Figure 2. Measured SDSS-II DR7 power spectra for red (top panel), blue (middle panel), and all (bottom panel) galaxies compared to 
several theoretical models. All curves plotted with thick lines use the cosmology from the WMAP5 ACDM+ ^2 m v MCMC chain that 
yields the maximum likelihood for WMAP5 only. Models | (iii) | and |(x)| are plotted with the bias parameters that give the best fit to the 
P (k) data shown for k < 0.3 /iMpc -1 and models [(ijj | (ii) | and | (vii) | are plotted with the best-fit b for the data shown at k < 0.1 /iMpc -1 
to better illustrate the different behaviours at small scales. The thin orange line and the thick dashed orange line both use model [(vii) | 
but the thin orange line uses cosmology and bias parameters that give the maximum likelihood to the P (k) data at k < 0.3 hMpc" 1 for 
the P ga i (k) = b 2 P p t (k) model - this cosmology has a high value of with a high and fits the data better than the thick dotted 

orange line (though not as well as models |(iii)| and |(x)| . This illustrates how increasing the neutrino mass can lead to a better fit for 
some models. See text for more detailed discussion. 



where the bias b is allowed to take on different values for 
different galaxy subsamples. 

(ii) Our second model replaces the linear matter 
power spectrum w ith the Halofit nonlinear transformation 
IjSmith et al.ll2003l ). which we denote Phf (k): 

P g ai (k) = b 2 P hi (k) . (6) 



(iii) For the next model we add one more free parameter 
to model [(!)] as an additive constant term: 

Pgai(fc) = b 2 P hn {k)+P . (7) 

The Po parameter has a physical basis: it represents a change 
in the shot noise. If galaxies are not a Poisson sample of 
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the matter distribution, our power spectrum calculation will 
subtract the shot noise imperfectly. Non-Pois son shot nois e 
is a generic prediction of the halo model (|Selialdl2OO0h : 
in this context, Po is an approximation to the one-halo 
contribu tion to the p o wer spectrum. This m o del was pro- 
posed bv lSeliakl (|200lD ; ISchulz fc White! <l2006l);lGuzik et all 



■ — « ■< ■ r- - n 

( 2007 ) and has been explore d by, e.g., IHamann et alj (J2008); 

ICresswell fc Percivall (|2009h . We refer to it as the P-model. 



(iv) We can do the same for model (i) giving 



P ga i(fc) = fc 2 P hf (fc) + P . 



(8) 



(v) As an altern ative to the P-mo del, we also test the Q- 
model put forth bv lCole et alj l|2005l ). which adds a different 
free parameter Q to "correct" for unknown scale-dependent 
galaxy bias and redshift space distortions: 



^gal (k) 



: — Pun (fc) . 



(9) 



1 + 1.4fc 

ICole et all (|2005h test and calibrate this model with N-body 
simulations populated by a semi-analytic galaxy formation 
model under the assumption of massless neutrinos - its va- 
lidity in the massive neutrino case has not been tested in 
detail. We include it here to facilitate comparison with othe r 
analyses that adopt this model, e.g.. iTegmark et all {2006). 

(vi) It is not clear that either the P-model or the Q-model 
correctly accounts for the damping of the Baryon Acoustic 
Oscillations (BAO) by the non-linear structure growth and 
redshift-space distortions. In order to confirm that this is not 
a big effect for neutrino mass constraints, we also consider 
a model with 



Pgal (fc) = fe2 TTrifc Pdw(fch 



(10) 



where Pj w (fc) is the "dewiggled" powe r spectrum calculated 
in the manner of lTegmark et all (|2006h that accounts for the 
non- linear damping of the BAO: 

Pdw (fc) = W (fc) Pita (fc) + [1 - W (fc)] Pnowigglc (fc) , (11) 

where W (fc) = e~' k//k *' ^ 2 and P n0 wi g gio (fc) represents a 
smooth power spectrum without BAO that we calculate us- 
ing a smooth b-spline fit to k Pun (fc) with eight nodes 
equally spaced in log (fc) from fc = 0. 007 to fc = 0.7 Mp c -1 . 
This is similar to the method used in Reid et al. (2010b) for 
massive neutrino models. Following ITegmark et al.l (j2006) 
we define ft, = a^^'^a^ (j4 s /0.6841) _1 ^ 2 with a± and cry 
given by equations (12) and (13) in lEisenstein et ail l|2007h . 

(vii) We now consider a model where we replace the linear 
matter power spectrum in (i) with P pt (fc) given by equa- 
tion ©: 



Pgal (fc) = 6 2 Ppt (fc) : 



(viii) We can do the same for the P-model (iii) 



Pgal (fc) = b 2 Ppt (fc)+P . 



(12) 



(13) 



(ix) For our next model, we use an extension of P pt (fc) to 
include nonlinear galaxy bias. Startin g from the assump tion 
of local deterministic nonlinear bias, ISaito et all (|2009l ) de- 
rive an expression for P ga i (fc) using the perturbation theory 



methods from iMcDonaldl (|2006l ): IJeong fc Komatsul (|2009l ) 
involving one new free parameter b 2 : 

P ga i (fc) = b 2 [P pt (fc) + b 2 P b2 ,s (fc) + b%P b22 (fc)] + Po, (14) 

where Pb2,s (fc) an d Pbii (fc) are given by equation (41) in 
ISaito et al. I (|2009h and b 2 is a third free parameter allowed 
to take on different values for different galaxy subsamples. 
We label this model as "PT with b, Po, b 2 " in our results. 

(x) Next we consider extending the pertu rbations the- 
ory a nalysis into redshift-space following iTaruva et alj 
(2009 ) and applyin g the redshift s pace distorti o n mod el of 
IScoccimarrol l|2004l) . According to IScoccimarrol (|2004l ). the 
redshift space power spectrum of matter can be modeled by 



Pzdist (fc, ju) = [Pss (fc) + 2f^ 2 Pse (fc) + / V Pee (fc)] 

/ ✓ - , v OA 



x exp (— (ffxka v ) 2 ) 



(15) 



where / = d\nD/d\na (a being the scale factor and D be- 
ing the growth factor) is the linear growth rate , fj, = k • z 
is the cosine of the angle between the wavevector fc and the 
line of sight z. Pss (fc) and Pee (fc) are nonlinear density and 
velocity power spectra, and Pse (fc) is the density-velocity 
cross-power spectrum, all calculated with one -loop pertur- 
bation theory as defined in equations (63-65) in lScoccimarrol 
l|2004l) . The one-dimensional velocity dispersion a v is defined 



1 f d 3 q Pee(g) 
3 7 (2tt) 3 q 2 ' 



(16) 



Following ISaito et all l|2009l ). we calculate Pss (fc), Pee (fc) 
and Pse (fc) in the same manner as in equation @: we ap- 
ply the one-loop corrections to the cold dark matter plus 
baryon portion of the power spectrum and leave the neu- 
trino power spectrum linear. Thus Pss, pt (fc) = Ppt (fc) as 
defined in equation (|3]), and Pee.pt (fc) and Pse, P t (fc) are de- 
fined analogously. 

To model the angle-averaged redshift space power spec- 
trum, which is what our data actually represents, we inte- 
grate equation (|15|) over ji: 



Pzdist (fc) = A) (fc) Ppt (fc) + 2fA 2 (fc) P M ,pt (fc) 
+f 2 A 4 (fc) Pee.pt (fc) 



(17) 



where A n (fc) are moments of the Gaussian term in equa- 
tion ©: 

A n (fc) = j dn^i" exp (- (/ nka v ) 2 ) (18) 

We model the galaxy power spectrum in redshift space with 
linear bias by using bias b for the density terms and no bias 
for the velocity terms, giving 



Pgai (fc) = A (fc) & 2 P pt (fc) + 2fA 2 (fc) WVpt (fc) 
+f 2 A 4 (fc) Pee.pt (fc) 



(19) 



This model is labeled as "PT with b, zdist". 

(xi) The P-model under this redshift space distortion 
model becomes 

Pgal (fc) = A (fc) [& 2 P pt (fc) + Po] 

+2fA 2 (fc) bP M ,pt (fc) + fAi (fc) Pee.pt (fc) (20) 
which we label as "PT with b, Po, zdist". 
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(xii) Finally, we combine the nonlinear bias model of 
equation (I14|l with the redshift space distortion model to 
give 

1.2 



Pgai (k) = A (k) {b 2 [P pt (k) + b 2 P M ,s (k) 

+b 2 2 P b22 (*)] +P } + 2A 2 (k) b [fPse.pt (k) 
+b 2 P b2v (k)] + f 2 Ai {k) Pee, pt (k) , 



(21) 



which we label as "PT with b, Po, b 2 , zdist". P b2v (k) is de- 
fined analogously to P b2 ,s (k) by replacing the density ker- 
nel wi th the velocity kernel in equation (41) in lSaito et al.l 
(l2009h . 

We calculate all of the various power spectra in equa- 
tions lfl"9)) . (f20)) . and (f2Tj) using a modified version of CAMB 
provided by Saito (private communication). Note that our 
redshift space distortion model does not add any new free 
parameters: P z dist (k) is fully defined by equation (|17[) once 
Pin (k) has been determined. 



Models (i 
ter b, models j(iii 



(ii) ( 



and (x) have one free bias parame- 
viii) | and | (xi) | have two free 



bias parameters {b and either Po or Q), and models (ix) and 
(xii) have three free bias parameters (6, Po and b 2 ). These 



free parameters are all allowed to take on different values 
for the red and blue galaxy subsamples. 

To illustrate some of the main features of these models, 
we plot a subset of them with our power spectrum data 
in Fig. [2] For all but one of the curves plotted, we use 
the parameter set from the ACDM+J^ m v MCMC chain 
that gives the maximum likelihood for WMAP5 only: l =( 
Qa = 0.7263, D. b /Q m = 0.1666, h = 0.7050, a 8 = 0.8065, 
n s = 0.9580, ^2m v = 0.0287 eV). The remaining curve is 
plotted with the parameter set from the MCMC chain giv- 
ing the maximum likelihood the all-galaxy P (k) data using 
the b 2 P pt (k) model: 1 = {Q A = 0.5168, D. b /Q m = 0.1392, 
h = 0.5596, cr 8 = 0.5459, n s = 0.9140, £ra„ = 1.34 eV). 

Figure[2]shows that model [(xj| (PT with 6, zdist) closely 
mimics model (iii) (the P-model) for the maximum likeli- 
hood WMAP5 cosmology and in fact fits the data just as 
well with one less free parameter. The other models plot- 
ted show different shapes at small scales, with the6 2 Pi n (k) 
model underestimating the small-scale power and b 2 P pt (k) 
overestimating it. The b 2 P pt (k) curve for the high ^ m v 
cosmology gives a shape closer to the data than the same 
model with WMAP5 cosmology because it compensates for 
the excess of small-scale power by adding more massive neu- 
trinos that damp it out. This means that a model of bias or 
nonlinearity that overestimates the small-scale power could 
potentially lead to a "false positive" measurement of neutino 
mass. The effect of this can be seen quantitatively in our re- 
sults. 



3.3 Calculating the galaxy power spectrum 
likelihood 

To compare our theoretical models to the power spectrum 
data, we convolve our model P ga i {k) with the window func- 
tion W {kj, k n ) discussed in Section [3.11 where kj are the 
central wavenumbers of the observed bandpowers and k n 
are the wavenumbe rs at which we evalu ate the theory power 
spectra. Following iReid et ail (|2010bT l. we define the con- 
volved theory power spectrum by 



Pgal,win {kj , 0, b ) — 

W {kj, k n ) P gal (k n /a scl , 0, b ) /aL - W {kj, 0122) 

71 

where are the cosmological parameters defined in Sec- 
tion [22] and b are the galaxy bias parameters for the model 
being used: b and possibly b 2 and Po or Q. 

Equation (|22p includes two very small adjustments: the 
W {kj , 0) term accounts for the integral constraint issue 
from e stimating the average galaxy density fr om the sample 
itself i|Percival et all 120071 : IReid et al.ll2010bh . and the a sc i 
factors adjust for the fact that we use one fiducial cosmology 
to convert our galaxy redshifts into distances rather than re- 
calculating the distance s for each cosmology (Tcgmar k et alj 
120061 : IReid et al.ll2010bl ). with 



D\ ( 2 mcd) 



Hi 



fid l^mcd) 



\D^{ Ziaed )?H{z mei ) 



1/3 



(23) 



where Da and H are the angular diameter distance and 
Hubble parameter for the cosmology defined by 0, D A d and 
.fffid are for the fiducial cosmology, and z me d is the median 
redshift of the galaxy sample. 

The likelihood was calculated assuming that the power 
spectrum data are distributed as a multi-variate Gaussian 
distribution: 

- 2\nL{0, b ) = x 2 {0, 0b) =X>jC7k A fe> ( 24 ) 

jk 

where Aj = [p ga i {kj) - P ga i,win {kj, 0, 0i,)j , P ga i {kj) is the 
measured galaxy power spectrum, and Cjk is the covariance 
matrix calculated as described in Section [3. II 

Our final galaxy power spectrum likelihood is defined 
by marginalizing over the bias parameters b by integrating 
the likelihood over the prior distribution n {0 b ): 



L 



p(k) 



{0) 



d n b L{0, b )n{0 b ) 



(25) 



We use a flat prior on b 2 from to +oo and flat priors 
on Po, b 2 , and Q ranging from — oo to +oo. We perform 
this marginalization analytically where possible, which is the 
case for most of our models. For models |(ix)| and |(xii) 
have to integrate over b 2 numerically, and for models 



(xi) and (xii) (the ones including redshift space distortion) 
we must do the integral over b 2 numerically as well. 

3.4 Importance sampling 

Because we have to evaluate the neutrino mass limit for 
several different power spectrum models and k cutoff values, 
we reduce the computational r equirement by using t he tech- 
nique of importance sampling (|Lewis & Bridle 2002) : we use 
a fixed WMAP MCMC chain and reweight the chain entries 
according to the revised likelihood values calculated for the 
different galaxy power spectrum fits. Given an MCMC chain 
of parameter values drawn from a likelihood distribution 
L, one can compute parameter constraints relative to a sim- 
ilar distribution L' by re-weighting the sample according to 
the likelihood ratios: 



L' (0*) 



-w L . 



(26) 
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Since the CMB measurement and the galaxy power spec- 
trum measurement are independent experiments, we can 
find the combined likelihood simply by multiplying the two 
likelihoods together: 



Lit 



CMB 



-p(fc) — Ap(*o ( l ) L\ 



0' 



Thus the WMAP MCMC chain can be reweighted by 



W CMB+P(fc) 



p(k) 



8' 



WCMB, 



(27) 



(28) 



so the 95% confidence upper limit on ^ m„ is the value of 
M" m that satisfies 



.95 = 



M*<M}>' 



J P{k) 



0' 



N 



T * 

1 P(k) 



(29) 



As a test of the importance sampling method, we applied 
it to the dark matter halo power spectrum Phaio{k) derived 
from th e DR7 Lum i nous R ed Galaxies using the likelihood 
code of iReid et all (|2010bT ) to evaluate Lpqa Our impor- 
tance sampling method yields Yl m v < 0.64 eV (95% CL) 
for WMAP5+Ph a io(fc), w hich agrees within 3% to the result 
from IReid et al.l l|2010br i of £ m„ < 0.62 eV. We therefore 
conclude that importance sampling is sufficiently robust for 
the purposes of our analysis in this paper. 

To summarize, the steps to calculate the neutrino mass 
limit from WMAP5 + galaxy power spectrum are: 

(i) For each entry of cosmological parameter values in 
the MCMC chain, calculate a linear power spectrum Pn n (fc) 
with those parameters using CAMB. 

(ii) Choose a model for nonlinear evolution and bias, as 
described in ij3.2l 

(iii) Calculate likelihood of this model power spectrum 
using a measured galaxy power spectrum. 

(iv) Marginalize over the bias parameters by integrating 
the resulting likelihood analyticall y or numerically . 

(v) Use importance sampling ( Lewis fc Bridle! 120021 ) to 



reweight the entry in the WMAP5 MCMC chain by the 
likelihood from the galaxy power spectrum: 



^CMB + P(fc) — ipfkl^CMB 



(30) 



(vi) Repeat this for each entry in WMAP5 MCMC chain, 
and then calculate the neutrino mass limit using histograms 
of the MCMC chain weighted by Wcmb+p(/c)- 



4 RESULTS 

4.1 Neutrino mass constraints 

We have performed the above procedure using the power 
spectra for the full SDSS-II DR7 main sample, and for the 
same sample split into red and blue galaxies. For each of 
the three samples, we have done the calculation using power 
spectrum values for k < fcmax, with fcmax ranging from 0.06 
to 0.3 h Mpc" 1 and repeated this for each of our twelve mod- 
els for bias and nonlinearities. Results for the first six models 
are plotted i n Fig . [3] and results for the six models based on 
ISaito et al l (|2009h are plotted in Fig. f4]The neutrino limits 
for the all-galaxy sample typically trace the limits for the 
red galaxies, indicating that the red galaxies dominate the 
signal - this is because the red galaxies are more highly bi- 
ased. Furthermore, the limits from the blue galaxy sample 



tend to be larger than those from the red galaxies — this is 
to be expected since the measurements using the red galax- 
ies are intrinsically more sensitive and therefore give more 
stringent limits. We also typically obtain weaker limits as 
we increase the number of free parameters used to model 
galaxy bias and nonlinearities, as expected. 

One notable feature seen in Figs. [3] and [4] is that for 
the all-galaxy sample using the 6 2 Phf (fc) model and all three 
galaxy samples using the b 2 P p t (fc) model, the neutrino limit 
increases sharply at the high end of our fc max range. To 
explore what is driving this, we plot the histograms for 
fcmax = 0.3 ft Mpc -1 for the red galaxy sample using the 
[PT w/ 6, zdist], b 2 P hf (fc) and b 2 P pt (fc) models in Fig. [5] 

The distribution for the b 2 P p t (fc) case clearly peaks at 
a non-zero neutrino mass, and a hint of a secondary peak 
at non-zero neutrino mass can be seen in the b 2 Phf (fc) his- 
togram as well. From Fig. [2] we can see that the driving force 
behind this effect is that for models with no flexibility in the 
power spectrum shape, the differences in the amount of up- 
turn at the high-fc end are compensated for by changing the 
cosmology instead. For the b 2 P pt (fc) model, the shape of the 
power spectrum for the best-fit WMAP5 parameters over- 
estimates the data at high fc, and increasing the neutrino 
mass can pull the high-fc end down into better agreement 
with the data. 

We have compared the neutrino mass limits calculated 
from the SDSS-II data against those from 200 mock cat- 
alogues: the shaded regions in Figs. [3] & [4] contain 68% 
of the neutrino mass limits from mock catalogues. As dis- 
cussed in Section 13.11 we would expect to measure the 
neutrino mass limit to be within this expected sensitivity 
range if the neutrino mass were equal to zero and we un- 
derstood bias perfectly. They will therefore underestimate 
the true range of limits. Because the mock catalogues are 
constructed to have a power spectrum matched to Py ln (fc) 
with linear bias and no redshift space distortions, we use 
the b 2 P\i n (fc) model on the mock catalogues to define the 

and (x) 



expected sensitivity range for models (i) (ii) (vii 



(all of the one-parameter bias models) . Likewise, we use the 
b 2 Pn n (fc) + Po model for the expected sensitivity range for 
and 



(viif 



(all of the P- models ), the 



models (iii) | (iv) 

b 2 (1 + QFj7 (! TTJk) PnJk) model for both of the Q- 
models (v) and (vi) and the [PT w/ b, Po, b%\ model for 
both of the 62 models (ix) and (xii) This correctly captures 



the weakening of sensitivity due to adding more free param- 
eters but still allows for a good fit to the power spectrum 
built into the mocks. (Note that this is not strictly true for 
the Q- models as Q = does not reduce to Pn n (fc) ; however, 
we expect the effects of this to be negligible.) 

Comparing the data curves to the shaded regions in 
Figs. [3] and [4] we see that the results from red and blue 
galaxies are generally consistent with the limits predicted 
from the mock catalogues at the la level, which suggests 
that current neutrino mass limits from large scale structure 
are sufficiently robust to differences in bias modeling. How- 
ever, the constraints from the red galaxies are systematically 
better than expected, while the results from the blue galaxies 
are systematically worse. We consider this difference further 
in Section \4. 31 

We show our 95% confidence limits on the summed neu- 
trino mass for each model Table [T] using both conserva- 
tive (fcmax < 0.12 /iMpc -1 ) and aggressive ranges of scale 
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Figure 3. 95% CL limits on ^2m v as determined by WMAP5 plus the SDSS-II DR7 g al axy power spectr u m up to a maximum 
wave n umber fc max for bias models frequently used in the literature, e.g. lEIgar0v et al.l j2002t) ; ICole et al.l d2005h ; iTegmark et al'. I l l2004 
2006); Hamann et al. | J2008h . Expected sensitivity ranges show the mean limit iW^, lm ± one standard deviation measured from 200 mock 
galaxy catalogs and represent what we would expect to measure for the neutrino mass limit assuming that m u = and we understood 
the bias perfectly. Green points indicate cases where the CMB and galaxy datasets may potentially be inconsistent with each other 
- these points should be regarded with caution. Using data at larger k values leads to tighter limits on Y2 m v but makes the results 
increasingly model-dependent. 
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Figure 4. 95% CL limits on ^2 m„ as determined by WMAP5 plus the SD SS-II DR7 ga l axy p o wer spectrum up to a maximum 
wavenumber fc ma x for bias models based on perturbation theory results from IScoccimarr 3 \2O04) ; ISaito et al.l <200Sft ; iTaruva et al.l 
(2009). Expected sensitivity ranges show the mean limit M 1 ™ ± one standard deviation measured from 200 mock galaxy catalogs and 
represent what we would expect to measure for the neutrino mass limit assuming that m„ = and we understood the bias perfectly. 
Green points indicate cases where the CMB and galaxy datasets may potentially be inconsistent with each other — these points should be 
regarded with caution. Using data at larger k values leads to tighter limits on ^2 m v but makes the results increasingly model-dependent. 
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Figure 5. Probability distributions for M v using WMAP5 + P(fc) with fc < 0.3 /iMpc 1 for the red galaxy samples. 

Table 1. Summary of neutrino limits, calculated to a conservative wavenumber limit fc max < 0.12 feMpc" 1 , and to a more aggressive 
limit fc max < 0.2 feMpc -1 . Models are classified by the number of free parameters that are used to model galaxy bias. 



Limits on ^ m v in eV at 95% CL 
Model 



^ of free Conservative (fcmax < 0.12 hMpc 1 ) Aggressive (fcmax < 0.2 fcMpc" 1 ) 





params 


blue 


red 


all 


blue 


red 


all 


& 2 Pli„ (fc) 


1 


0.98 


0.58 


0.57 


0.68 


0.45 


0.56 


b 2 P hf (fc) 


1 


1.00 


0.65 


0.65 


0.69 


0.51 


0.68 


b 2 P lin (fc) + P 


2 


1.32 


1.10 


1.03 


0.92 


0.73 


0.88 


6 2 Phf (fe) + Po 


2 


1.31 


1.10 


1.03 


0.92 


0.72 


0.85 


b 2 (l + Qfc 2 )/(l + 1.4fc)P lin (fc) 
fc 2 (l + Qfc 2 )/(l + 1.4fc)P dw (fc) 
6 2 P P t (fc) 


2 
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1.24 


0.97 


0.97 


PT with b, zdist 
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0.96 
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0.46 
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2 


1.31 


1.10 


1.03 


1.06 


0.75 


0.74 


PT with b, P ,b 2 , zdist 


3 


1.38 


1.17 


1.10 


1.25 


0.99 


0.97 



I 

iii) 

iv) 

Mi 

33, 



(fcmax < 0.2 ftMpc -1 ). The effect of the number of free pa- 
rameters in the model (excluding the effect of neutrinos) 
clearly has a strong effect on the recovered constraint. This 
is discussed further in Section [5] 



4.2 Consistency between CMB and LSS data 

When combining results from independent datasets such as 
our combination of CMB and galaxy power spectrum data, 
it is important to check that they are consistent with each 
other under the model. If the galaxy power spectrum is not 
fit well by any set of parameters that has a reasonably high 
likelihood with respect to the CMB d ata, results obtained by 
combining them will be meaningless. iMarshall et alj (2006) 
details a method to determine whether multiple datasets 
are consistent: compute the ratio of Bayesian evidences for 
two hypotheses Ho and Hi, where Ho is "There is one set 
of parameters defining one global model that describes all 
datasets," and Hi is "There is a different set of parameters 
that describes each dataset." 

For two independent datasets di and d,2, this ratio is 



(see equation (3) of IMarshall et all (|2006f l) 

Pr(di, d 2 \H ) 
Pr(di\Hi)Pr(d 2 \Hi) - 



R : 



(31) 



Because calculating the evidence requires integrating over 
the prior distribution for all of the parameters of the model, 
it is computationally difficult to calculate. Motivated by 
iKass fc Raftervl (| 19951 ). we define an approximation to the 
evidence that can be easily computed from an MCMC chain: 



Pr(di) 



(32) 



This is the harmonic mean of the likeliho od values in the 
MCM C chain and can be derived following IKass fc Raftervl 
(1995) by using importance sampling: since the MCMC 
chain samples the posterior distribution, we can approxi- 
mate an integral of Pr(di|0 8 ) over the prior distribution 
with a sum over the MCMC chain by reweighting using 
equation (|26|l with L' equal to the prior distribution 7r (0 1 ) 
and L equal to the post erior distribution Pr (fl'|rii) = 
Pr(di|0 i )7r(f? 1 )/Pr(d 1 ). l|Kass fc Raftervl 119951 ) indicate 
that while this quantity can be unstable in some cases, it 
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is typically accurate enough for making rough categoriza- 
tions, which is how we use it here. 

Following this same idea, we derive an approximation 
for R in the case where we have an MCMC chain for d\ and 
are using importance sampling to combine d\ and d 2 , as we 
describe in §j33J) for di=WMAP5 and d 2 = P ga i(fc) : 

R = z r f r , (33) 

/{ZtM pr(dil^)]- 1 } 

where w\ +2 is given, e.g., by equation (|28p . This is quite sim- 
ilar to the heuristic approach proposed by I Lewis fc Bridle! 
(2002) of comparing the mean likelihood of the samples 
where the posterior distributions overlap to the mean like- 
lihood under the original posterior: here we have sim- 
ply replaced the me an by the harmonic mean. Following 
iMarshall etakl (|2006l ). we interpret a value of R > 1 as indi- 
cating that the null hypothesis Ho of a global parameter set 
is favored over Hi with a separate parameter set for each 
data set. Conversely, R < 1 indicates that Hi is favored. Ob- 
taining a value of R < 1 does not automatically guarantee 
that the datasets disagree, but it is an indicator that a pos- 
sible mismatch between the datasets should be investigated 
further. 

In this analysis we use R as a conservative "warning la- 
bel" : we mark models where R < 1 as an indication that we 
should proceed with caution in interpreting the results since 
there may be disagreement between WMAP5 and P ga i(fc) 
under these models. These points are clearly marked in 
Figs. 0& 21 In particular, this warning label applies to the 
points at the high fc max end of the 6 2 P pt (k) results where 
the measured limits from the combined datasets are higher 
than the WMAP5 only limit. This suggests that applying 
such a consistency check to future datasets could avoid a 
spurious detection of neutrino mass caused by an incorrect 
bias model. 



4.3 Comparison between results from blue and 
red galaxies 

To display these effects more clearly, and include the cor- 
relations between the blue and red galaxy samples (caused 
by the overlapping volume covered), we plot the difference 
between the neutrino mass limits from red and blue galaxies 
as compared to the expected range from the mock catalogs 
in Fig. [6] To define the expected range, we compute the dif- 
ference between the neutrino limits derived from the mock 
blue and red galaxy samples drawn from the same mock LN 
density distribution. Then we find the mean and standard 
deviation of this difference over the 200 mocks, and plot the 
expected range as the mean ± one standard deviation. This 
accounts for correlations between the red and blue galaxy 
samples due to being drawn from the same underlying dis- 
tribution. 

Figure [6] shows that the difference between the blue and 
red galaxy limits is consistent at approximately the la level, 
compared with the expected range for nearly all of the dif- 
ferent bias models the entire fc max range tested. Red galaxy 
limits are systematically higher than those from the blue 
galaxies, even after accounting for the expected difference 
in sensitivity, but the significance of this is low. 



5 CONCLUSIONS 

We have performed a quantitative analysis considering 
how large-scale structure neutrino mass limits derived from 
galaxy surveys depend on our ability to model non-linear 
effects such as galaxy bias and redshift-space distortions. 
Galaxies form a bimodal distribution in colour and can be 
split into red and blue galaxies. These populations cluster 
in very different ways, with red galaxie s being more strongly 
clustered on large s c ales than b l ue (Norberg et al.l 120021; 



Tegmark et al.l 120041; I Wild et all 120051 ; ICollister fc Lahavl 
20051 ; IS wanson et al.l 120081 ), a nd have a different redshift - 
space power spectrum shape |Cresswell fc Percivall 120091 ) . 
Splitting into red and blue galaxies and analysing the sam- 
ples separately is therefore an obvious way to empirically 
analyse the robustness of neutrino mass limits from galaxy 
clustering measurements, and we have performed this test 
using data from the SDSS-II DR7. To evaluate the robust- 
ness of the derived neutrino mass limits, we compare the lim- 
its to the expected sensitivity range calculated from mock 
Log-Normal catalogues: given the accuracy with which we 
have measured the galaxy power spectrum, the expected sen- 
sitivity range defines the neutrino mass limit we expect to 
measure if the neutrino mass were equal to zero and we un- 
derstood bias perfectly. 

Galaxy bias and redshift-space effects alter the mea- 
sured galaxy power spectrum from that of the matter power 
spectrum and many different models have been suggested 
to account for this. To see how the choice of model affects 
the results for the neutrino mass limit, we have considered 
12 of these models, described in Section 13.21 that span the 
space of currently-used models and al so include new models 
l|Saito et al.ll2009l ; iTaruva et al.ll2009l ). 

Our models can be classified by the number of galaxy 
bias parameters used: the simplest models use only one pa- 
rameter, namely the constant linear bias b which simply 
renormalizes the power spectrum. More complicated mod- 
els use two or three parameters to account for more complex 
scale-dependent effects. It is clear from the results presented 
in Figs. 13 &U that most one-parameter models break down 
for fc max gl 0.1 — 0.2, giving results beyond these limits that 
are inconsistent with the expected sensitivity range. The 
exception is mo del (x) (a combin ation of the perturbation 
theory model of Saito et al. 20091 an d the redshift space dis- 



tortion model of [Scoccimarr o 2004]), which fits impressively 
well with only one free parameter. 

All models with two free parameters give remarkably 
similar results: clearly the inclusion of the extra parameter 
allows all of these models to mimic the broad features of 
galaxy bias and redshift-space distortions, leading to consis- 
tent constraints. The only 2-parameter model that stands 
out is the dewiggled Q-model, which is unable to fit the 
shape of the blue galaxy power spectrum. Three-parameter 
models give self-consistent but not very informative results: 
the inclusion of another free parameter significantly de- 
grades our ability to see the signal of massive neutrinos. 
We also note that adding more free parameters increases 
the possibility of degeneracy with cosmological parameters, 
which could pote ntially bias the results for large k values 
|Saito et alj|2010h . 

Overall, we find constraints from the red galaxies that 
are slightly stronger than expected: they are at the la edge 
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Figure 6. Difference between M]) m as determined by WMAP5 plus the SDSS-II DR7 red and blue galaxy power spectra. Shaded region 
shows the mean Af ^ lm (blue) — Mj? m (red) ± one standard deviation measured from 200 mock galaxy catalogs and represent what we 
would expect to measure assuming that M v = and we understood the bias perfectly. 



of expected constraints, while those for blue galaxies are 
slightly less restrictive than expected, again approximately 
at the edge of the la range expected. In order to assess 
the significance of this potential discrepancy we need to 
allow for correlations between the recovered red and blue 
galaxy power spectra as they cover the same volume of the 
Universe. We account for this using mock Log-Normal cat- 
alogues, drawing overlapping red and blue mock samples 
from the same density fields. Figure [6] shows that the blue 
galaxy limit minus the red galaxy limit is slightly larger 
than expected: the difference in the neutrino mass result 
between SDSS-II red and blue galaxies is at the la limit. 
This matches the results shown in Figs. [3] & 13] indicating 
that the effect of correlations between red and blue power 
spectra is small. 



We therefore see that, with current SDSS data, recov- 
ered neutrino mass constraints are broadly consistent for 
red and blue galaxies. With only these data we cannot tell 
if the intriguing differences between blue and red galaxies 
are caused by noise, or if there is an element lacking in all 
of the bias models we have tested. Our results indicate that 
current cosmological neutrino mass constraints are generally 
robust to astrophysical systematic effects, but that analyses 
of future surveys will need to treat these effects carefully 
in order to set convincingly tighter constraints. It is clear 
that being able split a galaxy catalogue into red and blue 
galaxies provides an important diagnostic test of potential 
systematic errors, when measuring neutrino masses. When 
designing future surveys, it will be important to allow for 
such tests for systematics. 
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